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RESPONSE OF A PIASTIC CIRCUIAR PIATE TO 


A DISTRIBUTED TIME-VAHIING LOADING 

By 

Deene J. Weidman 
ABSTRACT 

The bending response of a plate to dynamic loading has been of 
interest in plasticity theory for some time. Generally, many restrictive 
assumptions are made on both the material behavior and the loading applied 
to the plate. The plate material is considered rigid-plastic, and four 
basic yield criteria are discussed and evaluated, with the Tresca yield 
condition being selected from these four criteria as allowing the most 
generality of solution. This dissertation, however, removes the restric- 
tions made concerning the applied loading for the Tresca yield condition. 
The plate is considered to *>e simply supported, and therefore, only one 
point (or circle) of discontinuity is considered to exist across the 
plate radius. The location of this circle of discontinuity is not 
generally constant but varies with time. The movement of this "hinge 
circle" is the key to solution of plastic plate problems under general 
loadings that vary with radius and time. The differential equation that 
defines the motion of the hinge circle is derived, and solved exactly 
for seme cases. A short computer program is also presented that allows 
the numerical solution of the general non-linear differential equation 
for hinge circle location. 

Many previous authors have considered the radial loading distribu- 
tion to be either a partially (or fully) uniform radial distribution or 


else a "concentrated" load, and either statically or impulsively applied 
in time. A few papers allowing slightly more general loadings are avail- 
able, but in all previous work the authors make assumptions that require 
the radial location of the hinge circle to decrease with time. This 
assumption simplified the solution of the problem considerably, but it is 
not a valid assumption in general. In the present approach, the hinge 
circle is allowed to move as the loading dictates, and outward hinge 
motion is seen to occur. In some cases of general loading, the hinge 
circle does only move inward, but the actual rate of deflection could 
not have been calculated by previous analyses. Eventually, however, the 
hinge moves to the center of the plate, and the plate deforms conically 
until the hinge circle finally disappears, leaving a rigid, deformed 
plate. All of the quantities of interest (plate velocity, bending moments, 
stresses, etc.) are written in terms of the variable hinge circle location. 
The location of the hinge circle is then defined in the general case. 
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INTRODUCTION 


The prediction of the plastic response of a simply supported circular 
plate has been a problem of considerable interest in recent years. Much 
of the current knowledge in plastic-plate problems stems from original 
analyses instituted at Brown University in the early 1950's by H. G. 
Hopkins and W. Prager and others. In these early papers, the plastic- 
plate material was considered to follow the Tresca yield criterion and the 
loading was considered to be uniform over an annular or circular region 
on the surface of the plate and statically applied. The load-carrying 
capacity (or load at which the initial yielding occurs) was determined 

and the rate of deflection was found. Even though these earlier papers 

\ 

(as examples, see refs. 1 through 6) solved for only a few of the quan- 
tities of interest in plastic-plate theory, these papers presented a 
basic approach to the analysis of general plasticity problems. 

After the first few papers on plastic plates were published and the 
significance of these analyses was realized, many more extensive analyses 
were initiated and other effects were investigated. The influence of 
changing the yield condition from the Tresca yield condition to other 
forms (von Mises, etc.) was investigated (refs. 4 and 7) and relatively 
small differences in carrying capacity were noted due to the change of 
yield condition (for the various loadings considered). 

The influence of including the inplane forces and thus allowing for 
membrane action, has also been investigated (refs. 8 through 14 ) and 
this effect is shown to cause the deflection to be reduced from the 
deflection due to bending theory. This is a very difficult problem 
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involving the interaction of several yield variables (two bending moments 
and two inplane forces). Sometimes the interactions between moments and 
forces are neglected (see refs. 13 and 14 for current interaction 
theories). These inplane forces are present in plastic-plate problems 
as soon as deflection starts, and become increasingly important as deflec- 
tions increase. Initially, bending action strongly predominates, indica- 
ting that membrane forces start to have a significant effect only when 
the plate central deflection becomes greater than about twice the plate 
thickness. Thus, if the initiation of deflection is of primary interest, 
then the need for inclusion of membrane forces is greatly reduced. These 
membrane analyses (refs. 8-l4) included small computational approxima- 
tions in addition to the assumption of a uniformly distributed impulsive 
loading. Thus, a bending and membrane solution for a plate under a 
loading that varies with radius and time is still not available, and the 
results presented for uniformly distributed impulsive loadings may not 
be directly applicable to more general loading cases. 

The consideration of additional material properties such as strain 
hardening, viscosity, elastic deformation etc., has also been made 
(refs. 7, 9, and 15 through 17). In these references, consideration of 
any one of these properties allowed the theoretical analysis to predict 
the experimental results well. For example , reference 7 shows that the 
addition of viscosity effects could yield agreement of the theoretical 
with experimental fined, maximum deflections; however, the magnitude of 
the viscosity coefficient is unknown, and its value was selected for the 
specific case of a uniformly distributed impulsive loading. 
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Interaction problems have also been receiving attention more 
recently (see refs. 18 through 21). Interaction problems are defined 
as those problems in which more than two stress resultants are con- 
sidered necessary to define the yield surface. The analyses including 
inplane forces (refs. 8 through l4) are examples of this type of prob- 
lem. A three- or four-dimensional yield surface is needed, with each 
stress resultant allowed a range of values before it causes yielding by 
itself. Within this range of values, an interaction exists between all 
of the allowed stress resultants. Most often, the interaction is deter- 
mined to be either linear or quadratic in these stress resultants. The 
influence of shear forces (see ref. 21) is a particularly important 
interaction problem that has not been widely analyzed. For example, 
even though it i s a well-recognized fact that shear forces are of primary 
importance for concentrated load prcolems, bending theory has still been 
used to predict the loading- carrying capacity of plates under a concen- 
trated load. Also, the success of using a viscous- shearing model of 
failure for ballistic impact problems (ref. 22) indicates the need for 
the consideration of shearing forces. 

The investigation of other boundary conditions (such as clamped 
supports or elastic supports, refs. 23 through 25) and other geometries 
(such as annular plates, refs. 26 and 27) has also been initiated. An 
attempt to discuss all of the literature on the plastic theory of thin 
plates would not be appropriate here. An excellent group of bibliogra- 
phies and reviews can be found in references 28 through 32. Those 
papers directly discussed and referenced herein are thought to be 
representative of the best of current work in this field. 
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Two experimental papers must be mentioned here as essential to 
proper justification of the theoretical analyses. Analyses employing 
piecewis linear yield surfaces indicate that points (or circles) should 
exist across which a discontinuity is expected. These discontinuities 
would not be expected for smooth yield surfaces. Such a discontinuity 
has been observed and recorded for a beam (ref. 33) > yielding the 
definition of a hinge. More significantly, such a "hinge circle" has 
also been recorded for a thin plate (ref. 3*0* These experimental 
observations lend creditability to the use of piecewise linear yield 
surfaces, which have only been considered approximations in the past. 

In nearly all of the papers referred to above, the radial variation 
of the loading was considered to be either uniform over all or part of 
the plate surface or concentrated at the plate center. The analysis of 
a more general radial load distribution is needed, and this variation 
is allowed in this thesis. 

The solutions of the references centered around determining that 
unique single locate c.i ■’n these plates at which the hinge circle occurs. 
Initially, static leadings were considered, but in a classic paper by 
Hopkins and Prager (ref. 2) the influence of a time variation was first 
considered. In this reference, a uniform load was applied over the 
surface of the plate for a finite length of time and then removed. This 
analysis introduced the fact that the location of the hinge circle does 
not remain fixed, but after the loading is removed, the hinge circle 
shrinks to a single point at the center of the plate. This moving 
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boundary separates two different regimes (or regions of solution) and 
causes many complications. The difficulty of the moving hinge circle 
is also important in impulsive loading problems, and therefore the 
problem of defining its motion had to be solved. 

After the initial papers were presented, the investigation of 
impulsive loadings soon proceeded to dominate the literature (see, for 
example, refs. 5, 6, 8, 11, 18, 27, 35, and 36). A general radial varia- 
tion of the loading on the plate has been attempted in only one case 
(ref. 36) where a Gaussian radial load distribution was assumed to be 
impulsively applied. 

The solution of the response of a plastic plate to a general time 
variation of loading is a considerably more involved problem. The 
location of the hinge circle varies strongly if the loading changes 
with time. Very few references are available that allow time variation 
of loading (for example, refs. 2, 37, and 38)* Considering only 
uniformly loaded plates, with a time variation that quickly decays 
toward zero and never increases, Perzyna (ref. 37) states that the 
actual shape of the time variation has very little effect on the final 
deflections. However, in this reference the uniform load is maintained 
on the plate for a considerable length of time before rapid removal 
of the load, and thus it might be expected to show only small variations 
due to load removal. 

Sankaranarayanan (ref. 38), on the other hand, has shown that 
for spherical caps under a similar uniform loading the time variation 
of the loading greatly influenced the final deformation. Therefore, 
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a need is seen for the general determination of the motion of the hinge 
circle, euid, indeed, this unique movement of the hinge circle is the 
key to solution of involved dynamic plasticity problems in plates. 

Although this dissertation is concerned with the dynamic loading 
of simply supported, rigid-plastic plates, the basic problem that 
prompted the analysis is the ballistic impact problem. This problem 
is concerned with determination of the response of a plate when impacted 
by a projectile (or spray of projectiles) traveling at a high velocity. 
Retaining the integrity of the plate is of prime concern. At low impact 
velocities, the response of the plate is entirely elastic. As the 
projectile velocity increases, the response of the plate becomes pre- 
dominantly plastic, and it is this problem of plasticity that is of 
interest herein. An excellent review of this impact problem (from the 
elastic point of view) is found in reference 39 > The basic point of 
interest is that the loading on the plate is definitely not impulsive, 
or indeed, not even uniformly distributed across the surface of the 
impacted plate. Therefore, a general analysis allowing a loading 


variation in both radius and time is desired. 



LIST OF SYMBOLS 


a,a 

b 

f(®i> 

f x (r,p,t) 


h 

i 1 (p,r),I 2 (p,r) 

M 

Mo 

P(t) 

Pp 

Pi 


Pmax 

Pi 

q(r) 

Q 


r,e,z 


R(r),T(t) 


parameters in example leading cases 
puter radius pf plate 

arbitrary functiens and cpnstants used fer selutiens 
functienal fprm pf the yield surface 
functicn used in scluticn fer regime AB, see 
equatien (50) 
half -thickness ef plate 

general functiens defined by equations (57) 
and ( 58 ) 

bending -mement resultant 

p 

yield -mement resultant, cr 0 h 
time variation of applied loading 
the value of loading at which initial hinge forms 
the value of loading at which hinge occurs away from 
origin 

maximum value of loading applied (maximum value of 
P(t)) 

specific points on a general yield surface 
radial variation of applied loading 
shear force resultant 

radial, circumferential, and transverse coordinates, 
respectively 

functions for separation of variables solution, see 
equation (4l) 
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s 

t 

t o 

t l 

t* 

t(r) 

t(p) 

u 

Vo 

V 

3 

€,€ 


4 

P(t) 

Po 

p l ,p 2 

a 

a o 

CT z ,T rz 

T o 

T r 0 ,T 0 z 


coordinate along p(t) curve, see figure 5 
time 

initial time 

time hinge circle disappears 

time plate comes to rest 

time needed for hinge to reach point r 

time needed for hinge to reach current location p 

radial displacement in plate 

arbitrary applied velocity for impulsive loading 
plate deflection 

general flow rule constant, see equation ( 33 ) 

strain and strain rate in plate, respectively 

dummy variables of integration 

curvature rate 

mass per unit area of plate 

general hinge circle location 

initial hinge circle location 

two general functions of p, see equation ( 65 ) 

an average stress, defined by equation ( 8 ) 

yield stress in simple tension or compression 

principal stress components 

vertical and shearing stresses, assumed small 

octahedral shear stress 

shearing stresses 
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Subscripts 

1,2 

r,e 


Extra Symbols 
[] 


refer to two planar coordinates (often r and 0) 
refer to radial and circumferential quantities, 
respectively 

denotes the jump in a quantity across the hinge 
denotes quantity evaluated at a given point 


GMERAL PLASTICITY CONSIDERATIONS 


In this section some of the basic considerations for any flow (or 
incremental) theory of plasticity are discussed. Several possible 
choices are evaluated, and from these choices the methods and cases con- 
sidered in the body of the text are selected. These methods have been 
chosen to yield a theoretical solution to a general group of plasticity 
problems . This is done without recourse to a large digital computer 
prograjn solving high-order, coupled differential equations and utilizing 
multivariable difference techniques, iteration procedures, etc., for a 
solution. 


Yield Criteria and Flow Rules 

Initially, a basic criterion of material yielding must be selected. 
The engineering material to be used, and the application to which it will 
ultimately be put, dictate the type of material yielding to consider. 

This philosophy of yielding (whether maximum tensile stress, maximum 
shear stress, etc.) defines the actual shape of the yield surface. Once 
the shape of the yield surface is determined, the generalized flow rule 
is applied. It is assumed herein that the material is always Isotropic 
during flow, so that principal stress directions are the same as the 
principal strain-rate directions . The flow rule then states that the 
principal strain rates are in the same proportions as the direction 
cosines of the outward facing normal to the yield surface at that point. 
This approach breaks down only in cases where the yield surface has sharp 
corners, and the outward normal direction is not unique. Such cases are 
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important, since many yield surfaces do possess corners. However, 
additional information is available in such cases that allows a unique 
solution to be found. This fact will be illustrated later in this 
section. 

In the case of an isotropic, perfectly plastic solid, the yield 
condition f(a i ) =0 can be considered as a three-dimensional surface 
in the principal stress space as shown in figure 1. The type of yield 
criterion selected would affect the yield condition f(a^) and thus 
change the shape of this surface. The flow rule states that a relation- 
ship exists between a given point (a^) on the yield surface and the 
principal strain rates (e^ ) at that point. The flow vector at a point on 
the yield surface has direction cosines with the ratios of values, 

and is normal to f(a^) = 0, or 


€ J " da/ 


( 1 ) 


where the quantity A is the constant of proportionality. 

The flow rule is then a function of the local normal to the yield 
surface and at any point on the yield surface at which the normal is 
unique (say points P^^ and P 2 in fig. l) the strain rate directions 
are well defined. Using strain displacement relationships or (taking 
time derivatives of these expressions) strain rate displacement rate 
expressions, the strain rates can be used to solve for the deflections 
in regions where the flow rule vector is well defined. However, at any 
point where the normal is not unique (say the point Pj in fig. l), the 
flow rule vector may lie anywhere within the cone of vectors that is 
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formed by the normals to all possible surfaces of f(a^) - 0 that pass 
through the given point. For the point P 3 , it appears that one flow 
mechanism corresponds to the surface containing , and another flow 
mechanism corresponds to the surface containing Pgj so that at the point 
Pj, a linear combination of these two mechanisms might be expected to 
occur. Thus, 3 times one flow mechanism is added to (l - 3) times the 
other flow mechanism, to yield the combined mechanism at the point Py 
The constant, 3, should then lie between 0 and 1. Note that the discon- 
tinuity between the P-j_ and P 2 types of flow mechanisms disappears 
toward the rearward portion of the yield surface, indicating that, in 
general cases, the occurrence of sharp corners may vary with the state of 
stress (or even with time). 

A short comment must be made concerning the influence of boundary 
conditions in yield criterion (or yield curve) selection. If a plane 
stress problem for a circular plate is assumed, the yield surface becomes 
a curve in the ct^ , a 2 plane. A portion of this plane is shown in 
figure 2. The three points (marked +) must be on the yield surface, 
and if a simply-supported plate is being analyzed, only the shaded region 
A is to be considered (otj. = Oq at the center of the plate, and a r = 0 
at the outer edge) . Any yield condition that has simple differential 
equations in this region is then of interest for simply-supported plates. 
However, if a clamped plate is being analyzed, both of the shaded regions 
A and B are of interest. An entirely different yield condition may 

have simpler differential equations in this larger region. 

In summary, the procedure to follow is to first select a yielding 
criterion and from this criterion determine the yield surface. Next the 
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specific associated flow rules are written for that particular surface, 
and the necessary equations for solution are written. In all cases, care 
must be taken to include the desired material properties, and even to 
consider the boundary conditions of a specific problem when selecting 
the yield criterion. Current yield criteria are discussed in the 
following sections to illustrate the manner in which this approach is 
applied. In these sections it is assumed for simplicity of discussion 
that the plate under analysis is in a state of plane stress ( = 0 ), 
and the yield surfaces are drawn as curves in the (a-pCg) plane. It is 
important to note that, while the yield surface is considered a curve in 
the oj = 0 plane, the normals to the three-dimensional surface are not 
in this (0^,02) plane. However, the projection of these normal components 
onto the (0^,0 g) plf'ne .is normal to the yield curve, and the principle of 
the flow vector is still applicable. Only bending stresses are considered 
herein, so that and are related, as well as op and • Thus, 

equations can be written in terms of either moments or stresses . Also 
it is a general fact that for an isotropic plastic body the yield surface 
is symmetric with respect to a ^ 5 ° line through the origin. With these 
general thoughts in mind, the specific yield curves will now be discussed, 
considering the x-^ coordinates are r, 9 , and z, respectively. 

Maximum octahedral shear stress - von Mlses .- The octahedral shear 
stress is written in general as 

x 0 = j\J(a 1 - o 2 ) + (a-L - 0j) + (a 2 - 0j) 


Ik 


and, for the plane stress case of a circular plate of thickness 2h 


„ _ M r 


m 0 

a 2 “ a e = 


and therefore 


a ° = (/f) T ° = (t§)\A* - a e> 2 + 4 + °e 

where a 0 is the yield stress in simple tension, and is related to the 
M 

yield moment °o = terms of the moment resultants, equation (2) 

can be written in the form of a yield surface equation as 

Mp + Mg - MyMg - = 0 (3) 

If equation (3) is plotted in the Mj., Mq plane, an ellipse (with its 
major axis at ^5°) is found as s.iown in figure 3* 

In the equilibrium equations developed later, the differential 
equation for Mj. contains a single term in Mq, and this term is 
removed by solving equation (3) for Mq and substituting into the 
equilibrium equation. From equation (3), then 

"e ' T i \l^ lh) 

The upper sign refers to the upper half (above the major axis) of the 
yield surface and the lower sign to the lower halt'. From equation (4) 
it is seen that the expression for Mq would cause a strongly nonlinear 
differential equation for Mj. with a nonremovable square-root term. 
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Therefore, any approach utilizing the von Mises yield condition would 
usually require a large-size computer program for solution of all the 
variables involved. However, it may be noted that this yield surface is 
the only yield surface considered herein that is smooth and regular (no 
corners), and various plasticity regimes (particular regions of solution) 
need not be defined for sue!' a yield function. 

Maximum tensile stress - Johansen .- In the case of piecewise linear 
yield surfaces, the definition of the yield surface (or yield function) 
is not a continuous function as it was for the von Mises yield surface. 
Instead, piecewise linearity means that the statement of the conditions 
for yielding should be expressed as a maximum or set of maximum state- 
ments. For a homogeneous material with a low yield strength in tension 
(such as a uniformly mixed concrete, for example), the tensile stress 
governs the yielding of the material, and the yield condition can be 
written In the form 


max. (1^1) - cr 0 = 0 (5) 

For the plane stress conditions considered herein, only the stresses 
a r and rr^ are included.. Then, if |cr r | is greater than | CTq j , then 
|°r! = 1^1 is the yield condition. This portion of the yield surface 
consists of two straight vertical lines at Oj. = ±a 0 , and M r = ±Mq on 
these lines. The differential equation for M r then just defines 
exactly Mg, and the solution is completed. Similarly, for |<jg| 
greater than | o r | , the yield surface becomes two straight horizontal 
lines at cr fl « ±a 0 and Mg ■ IMq along these lines. The differential 
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equation for Mj. is then solved exactly and the solution is again 
completed* The yield condition is therefore a square, and is shown in 
figure 4. The use of this yield criterion leads to simple, solvable 
differential equations to determine the moments and deflections. However, 
most metals have high tensile strengths and often fail in a shearing 
manner. Since the equations for the Johansen yield condition are quite 
similar in form to the equations for the Tresca yield condition (which 
is a shearing failure criterion) consideration of the Johansen yield 
surface is deferred here for future study. 

Maximum shearing stress - Tresca . - The maximum shearing stress yield 

criterion of Tresca is by far the most often used yield criterion. It 

defines a piecewise linear yield surface, and, therefore, the equation of 

the surface is expressed as a maximum type of statement. All points on 

Cf C 

the yield surface have a maximum shear stress equal to — , so that the 
equation of the yield surface (a six-sided surface in general) can be 
written 


max. 





J i i (6) 


Consideration of the plane stress problem analyzed herein (a^ = 0), and 
the yield condition of equation (6), the yield surface (in the a^, 0 % 

plane) becomes a hexagon as shown in figure 4. Since = a z = 0, the 
equations of the yield curve are directly written in terms of the moments 
as 
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h 2 a 0 = ^ = +M 0 h 2 a r = Mj. = *M Q My - Mg = *M Q (7) 

The upper signs refer again to the upper half of the yield curve, fech 
of these equations, by itself, can be easily solved for Mq, and this 
value of Mq substituted into the differential equation for Mj. . The 
equations for Mp are then solved, and the procedure for solution can 
be continued. However, each straight-line portion of the yield curve has 
its own unique solution, and it must be determined exactly which regimes 
apply to each part of the plate at each instant of time. Four of the 
corners of the Tresca yield hexagon would also be plasticity regimes of 
finite width on the plate. Thus, although the differential equations 
are simplified, the solution of any specific plate problem involves 
keeping track, in time, of each boundary between different regimes of the 
solution. Even with this difficulty it appears that use of the Tresca yield 
conation offers a consistent, well-known approach for the general solution 
of plate problems in plasticity, and therefore it has been selected as 
the yield condition to use in the main analysis of this thesis. 

Maxi m um reduced stress - Haythornthwaite .- This yield criterion was 
introduced some time ago (see ref. ^0), but has not been extensively 
investigated as yet. Since this yield condition is a piecewise linear 
yield surface, a maximum statement for yielding is expected, and the 
yield surface can be written as 


max. {[a i -a|}-|a o = 0 


(8) 
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(an + CTp + O5) 

where a is the average stress, or a = — . This six-sided 

surface is somewhat similar to the Tresca yield surface (see eq . ( 6 )). 

For the case of plane stress (05 = 0), equation ( 8 ) becomes a hexagon in 
the o±, a 2 plane as shown in figure 3 , and the equations for the sides 
become 


(2oi - ff 2 ) = *2a 0 (2a 2 - O].) = ±2<j q + a 2 = ±2a Q (9) 

where again the upper signs refer generally to the upper half of the 
yield curve, and the lover signs to the lower half. Since these stresses 
are related directly to the moments, equation ( 9 ) is seen to yield linear 
relationships between the moments Mj* and Mq as in the case of Tresca 
yield criterion. However, from figure 3 > it is seen that the Haythorn- 
thwaite yield condition closely circumscribes the von Mises ellipse, and 
can be used very effectively to provide, for example, an upper bound for 
load-carrying capacities of circular plates, while still retaining linear 
equations as in the Tr-sea yield condition. Also the Tresca yield con- 
dition can still be used effectively for a lower bound, and close approxi- 
mation to actual load -carrying capacities can be obtained. 

This yield condition also has additional advantages over the Tresca 
yield criterion in a few cases. The Tresca condition allows as many as 
three hinge circles to form, (when considering the upper half of yield 
surface), but the Haythornthwaite condition allows only one hinge circle 
to form. These differences can cause a considerable variation in the 
computational effort required for a general solution. In the analysis 
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of a simply supported plate, however, the region of interest (shaded 
region A in fig. 2) is such that the Tresca yield criterion allows 
either one or two solution regimes, whereas the Haythornthwaite yield 
criterion always requires two regimes. For this reason, the Tresca 
yield condition is assumed in the main body of this dissertation. 

Hinge Circles and Discontinuity Conditions 
The term "hinge circle" and the definition of such a cix-cl'-' for 
plastic-plate problems can be directly traced to beam analyses that were 
conducted previously. Experimental results for beams (ref. 33) have 
verified that points occur along the beam across which a large change 
in the character of the deflection exists. The location of these points 
varies with time, and propagates along the beam in dynamic loading 
problems. In fact, the measured velocity of propagation of these 
"discontinuities" along the beam can be compared with the expected rate 
of propagation for the "hinges" that naturally occur in the rigid -plastic 
analysis of a beam. The calculated velocity of the hinge (using the 
Tresca yield criterion) agrees quite well with the experimental propaga- 
tion rate except for very early times, and indicates that the rigid- 
plastic analysis yields reasonable results for a beam (see ref. 33)- 
Considering a generalization of this beam hinge concept for 
dynamically loaded circular plates, the term "hinge circle" was introduced 

to define the locations across the radius of the plate at which discon- 

1 ^ /^w\ 

tinuities in the circumferential curvature rate - p occur. These 

discontinuities have all of the appearances of hinges as in beams, and 
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again, experimental results (ref. 34) for plates generally verify the 
hinge propagation rates, even though the rigid-plastic theory yields 
infinite initial velocities under certain conditions. This agreement 
between experiment and rigid-plastic theory is somewhat surprising, since 
the material of the plate is inherently elastic, strain hardening, and 
under inplane and shear forces, etc., whereas the theory is derived for 
only the bending of a rigid-plastic material. This agreement in hinge 
circle motion implies that the bending theory gives a good estimate of 
plate motion until large deflections occur. 

Since moving circles of discontinuity are a fundamental part of 
dynamic plastic-plate pro lems, a definition is needed for the changes 
in the solution as a discontinuity circle passes any given point. These 
changes are the so-called "jump conditions" of plastic-plate theory, and 
can be determined from the knowledge of the quantities which are contin- 


uous . For example, to retain integrity of the plate, the deflection w, 


as well as the plate velocity — , must be continuous everywhere in the 

at 

plate, or 



( 10 ) 


where the bracket denotes the difference in the variable enclosed in the 
bracket between one side of the circle of discontinuity and the other 
side. Thus, the bracket denotes the jump of its argument as a discon- 
tinuity passes. Now, if the location of the circle of discontinuity is 
defined as p(t), this location can be plotted as in figure 5 as the 
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curve r = p(t) in a rectangular r, t space. Moving along this 
curve, the rate of change of deflection may be written 


dw _ dw dr | dw dt 
ds ds ds 


(ID 


Using the bracket notation presented in equation (in), the change of the 
quantities in equation (11) as the curve is crossed can be written 


' — — ' 
dw 


[dw 

§L + 

dw 

ds 


[dr^ 

ds 

^dt 


( 12 ) 


Since w is continuous along the curve, the left-hand side of 

equation (12) must be zero, and multiplication by ff- yields 

dt 


— ► — 

<5v 



dp(t) 


■ 

Sr 

w 

dt 


(13) 


This equation shows the basic form of the discontinuity jump conditions 


for plastic-plate problems. Since 


is also continuous everywhere, 

dw 


equation (13) simplifies considerably. Also replacing w by ^ in 
equation (13) yields another discontinuity condition; these two conditions 
can be written as 


Sfi 

dt 



= 0 


( 1*0 


"d£w d£ d 2 w _ n 

dt 2j + dt|d r dtj " 0 


( 15 ) 
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These two discontinuity conditions are consequei^?s of maintaining the 
continuity cf the plate (no breaking or plate failure). Dynamic equili- 
brium also requires that the moment M r and the shear force Q r be 
continuous in moving across a circle of discontinuity. By following the 
same procedure as for the deflections, the other two discontinuity 
conditions become: 


*M r 

+ 

dM r 

bt ^ 

dt 

Sr _ 



(16) 


In addition, for the case of a moving circle of discontinuity only, the 
only other second derivative jump condition can be written from equa- 
tion (l4). Since — is nonzero, ^ must be continuous, and therefore, 
crc a r 

from equation (13) directly 


f>wl 

+ &■ 


[dr dtj 

dt 

[dr 2 J 


(i * °) < i7 > 


With all four of the basic jump conditions derived in general, in 

equations (14) through (l6), the difference between moving and stationary 

discontinuities, and between hinge circles and nonhinge circles, becomes 

clear. If the discontinuity is not moving, then ^ is zero, and all 

dt 

terms not containing ^ must be continuous. Also, all bracketed terms 

dt 

that multiply * terms may then have discontinuities at the circle. 

U. X 
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If the circle of discontinuity is considered a hinge circle, then the 


terms containing 


are nonzero, and the other terms that appear in 


dr dt 

conjunction with r — ■*— must be discontinuous across the hinge circle, 
dr ot 

In general, whenever a bracketed term can be shown to be zero across 
some circle of discontinuity, then its argument may be substituted into 
equation (13) tor w, and another Jump condition on higher derivatives 
can be written. The basic Jump conditions, however, are still equa- 
tions (14) through ( 16 ). With these general jump conditions then, the 
solutions in two different regimes may be related, and problems solved 
more readily. If a problem to be solved contained two circles of 
discontinuity, the jump conditions that would apply, if the two circles 
ever met, would have to be derived for the case of intersection of 
two curves of the same type as the single curve shown in figure 5* After 
two such curves intersect, a new solution must be initiated, and the 
resultant single curve that would apply after intersection would 
necessarily have a discontinuity in its slope at the intersection point. 
Such a discontinuity in velocity of propagation would be bothersome to 
analyze, but since this dissertation concerns only simply supported plates, 
only one circle of discontinuity is possible, and therefore, the diffi- 
culties caused by interfering circles are not present. 


Plasticity Regimes 

After selection of the yield criterion and the boundary conditions 
for any given problem, the yield curve and its associated flow rule must 
be applied for solution of the problem. Depending on the selections 
made, several possible regions (or "regimes") of solution can occur. 
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The solutions in each of these regimes must be matched at the circle of 
discontinuity between them, so that the continuity (and discontinuity) 
conditions of the previous section are satisfied. In this manner, the 
complete solution is found. 

In the case of the von Mises yield criterion (see fig. 3) , and for 
either simply supported or clamped boundary conditions, the plate is 
always in a single regime with uniquely defined normal directions. 

Therefore, the problems of different regimes of plasticity and solution 
matching at a circle of discontinuity are avoided. However, the nonlinear 
equations are generally not solvable in closed form, and numerical 
integration methods must be employed for solution. 

All other yield curves discussed herein are piecewise linear, and 
various regimes of solution must be determined and matched. An additional 
difficulty also appears when the loading varies with time and radius . In 
very general cases, some regimes may appear and disappear as time increases. 
Therefore, proceeding to solutions in the general cases must be done 
slowly and carefully lest some regime of solution be improperly allowed 
or restricted. For piecewise linear yield curves, the linear portions of 
the curves are always possible regimes on the plate in the final solution. 
However, the corners of the yield curve may also be regimes of finite 
size, and each problem solved must be thoroughly analyzed for this 
possibility. Generally, it appears that a corner of the yield surface 
could represent a finite size regime on the plate only if that corner 
could possibly be a hinge circle. Indeed, if this is the case, then an 
estimate of the possible regimes for a given problem can be found simply 
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by the definition of a hinge circle as follows. Across a hinge circle, 

cV-ty 

- — — is discontinuous by definition, and is directly related to the 
strain rate. In fact, at the outer fibers of the plate 


dr St 



(18) 


If equation (1) is used at a hinge circle, then the normal to the yield 
curve has the direction cosines (e Q ) ; (e„). For unbounded e Q , the 
normal to the yield curve must therefore have the possibility of being 
horizontal. Thus, all corners of the piecewise linear yield curves that 
admit the possibility of having a horizontal normal are possible finite 
regimes on the plate. From figures J and b, for both simply supported 
and clamped plates, the possible finite corner regimes can then be read 
off directly: for Tresca, regimes A and D; for Johansen, regimes A 

and Cf and for Haythornthwaite, none. Table I shows a listing of all 
the possible plasticity regimes for transversely loaded, circular plates 
considering all four yield criteria discussed herein. Whether all (or 
indeed any) of these regimes occur in any particular problem depends upon 
the manner and magnitude of the transverse loading, and care must be 
exercised in the process of solution. 

Low Loading Cases - p < p < p n 

o - max. - *l 

As a general loading of sufficiently small magnitude is applied to 
a circular plate, the state of stress everywhere within the plate 
corresponds to a point inside of the yield curve, and for a rigid-plastic 


26 


material the plate must remain rigid. Irrespective of the type of loading 
and manner of application, the plate does not move or yield. As the 
magnitude of the loading increases, the plate has some states of stress 
corresponding to points closer to the yield curve, and, consequently, 
closer to yielding. As the magnitude of the loading increases still 
further, the state of stress at some location in the plate finally 
reaches a point on the yield curve and yielding starts . For a rigid- 
plastic material the stress is a a everywhere through the depth at that 
location, and the value of the loading that causes this yielding is 
defined as the "load -carrying capacity" of the plate. Since work 
hardening is neglected herein, all states of stress may not move beyond 
the original yield surface. The type of yielding that occurs and where 
it occurs is in the general case a function of the type of yield surface 
and the type of loading. In order to proceed further in this scussion, 
the Tresca yield condition for a plate with simply-supported edges is now 
assumed . 

The load at which yielding starts and a hinge circle is first formed 
is defined as p 0 . This value of loading is the load-carrying capacity 
for the Tresca yield criterion, which has been shown to yield a lower 
bound (ref. 4o) for the actual load -carrying capacity of a plate with 
any general yield criterion. Thus, p Q can be used as a conservative 
estimate for design of rigid-plastic circular plates. For the Tresca 
yield criterion, the value of p 0 is determined as follows: For initial 

yielding, a hinge circle (actually a point hinge) is formed at the origin 
of the circular plate. Then, Mj. = Mq = Mq, and the center of the plate 
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io located at point A of the Tresca yield curve shown in figure 4. 

At the edge of the plate, the moment Mj, is zero and point B of the 
yield curve applies. Therefore, the entire plate is in regime AB, the 
moment Mg * Mq everywhere in the plate, and (neglecting the inertia 
term since the rigid regime is just ending) the differential equation for 
Mj. becomes (see analysis section, eq. (39) )> 

^-(rMj.) - Mq = - J q(r)p(t)r dr (19) 


If the time axis is shifted so that t ■ 0 corresponds to the initiation 
of motion, then p(o) = p c . Solving equation (19) for Mj. then yields 
(using the condition that M r * M 0 at r = 0) 


Mr 


Mo • r Po 



q(r)r dr > drj 


} 


( 20 ) 


Nov, the boundary condition (that My at r = b is zero) is applied, 
and p Q can be calculated as 


Pc = 



( 21 ) 


This expression allows the calculation of the magnitude of p(t), for 
which a general q(r) loading will Just cause a hinge to form. If the 
loading is considered to be uniform over either a circular region around 
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the origin or an annulus near the supports, the expressions for load- 
carrying capacity shown in reference 1 for such loadings are found. 

For the problems analyzed herein, another quantity is also of 
interest. This quantity, called p^, is defined as the magnitude of the 
time variation of loading (p(t)) above which the hinge circle on the 
plate must « xist away from the origin. The determination of p^ is a 
much more difficult problem than the determination of p Q , and it has not 
been calculated previously except for an impulsively applied uniform 
load (ref. 2). If the loading applied to a plate is increased above the 
value p Q , then inertia terms become important in the differential 
equation as the plate continues to deflect. However, the hinge remains 
at the center of the plate, and therefore the single regime AB still 
applies everywhere on the plate. As the load is increased, 
eventually the load p^ is reached at which the distribution of moment 
Mp forces the hinge to move away from the origin. That is, below the 
load , the moment Mp = at the origin and decreases monotonically 
with radius to zero at the supports. Above the load p^ , however, the 
moment just beyond the origin attempts to reach a value greater than 
Mq . Since the yield criterion precludes such a value, the hinge, there- 
fore, must occur away from the origin. 

The value of p^ can be determined by solution of a general 
plasticity problem as follows: since the inertia force is important, 

the velocity must be determined In regime AB such that the velocity 
of the supports is zero-, anus, in regime AB, the curvature rate and 
support velocity are 
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K r = 




(22) 


Integrating twice with respect to r and using the boundary velocity 
equal to zero, the velocity can be calculated from equations (22) as 


§* = Ci(t)(r - b) 


(23) 


The differential equation for the moment M r in regime AB (with 
Mq - M Q ) is written (from the analysis section, eq . (39)) as 


^(i-Mr) 


Mo 


--/'f 

J o L 


q(r)p(t) 



r dr 


(24) 


Substituting the plate velocity from equation (23) into equation (24), 
and solving, the moment ILp becomes 


M - M EM 

M r - Mq - r 



(t - b)g d| 



(25) 

Two boundary conditions must be satisfied on Mj,. These are 

Mr. I = Mq and Mr I =0 (26) 

I r=0 1 r=b 

To apply the boundary condition at r = 0, the two integral terms in 
equation (25) must be evaluated using L'Hospital's rule and Leibniz' 
rule for differentiation under an integral sign. Thus, for example, 
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lim 
r -» 0 



lim 
r -> 0 



q(l)l d| 


(i) 


= 0 


(27) 


In this manner, both of the integral terms vanish, and application of the 
first boundary condition yields Cj = 0. Application of the second 
boundary condition is direct, and yields 



-bi^ + p(t) 



qU)l dg ) drj 



(28) 


Therefore, M r is fully known. Calculation of the derivative of M r 
follows directly from equations ( 25 ) and ( 28 ), and is 



By simultaneous application of L' Hospital ' s rule and Leibniz' rule, 
all of the terms in equation ( 29 ) become zero as r approaches zero. 
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Therefore, the origin is an extremum point for the moment. To determine 
whether the origin is actually a maximum point (as required by the yield 
criterion) or not, the second derivative of the moment must be evaluated 
at the origin. Again, the same method is applied, and at the origin the 
second derivative becomes 



(30) 


This expression changes sign from minus to plus at a value of p(t) 
equal to p-^, and equating expression (30) to zero yields 


Pi = 


bMo 

J q(i)i(b - i)d| - q -°^ 

d o 12 


(3D 


For values of loading less than this value, the moment curve has a true 
maximum at the origin. For values of loading greater than p-^, this 
moment distribution would increase with r from its value of M Q at 
the origin. Such an increase would violate the yield criterion and 
cannot be allowed. In such cases, the hinge circle moves outward from 
the origin, causing a finite region of regime A to appear, and the 
difficulties of a moving hinge circle must be surmounted. 
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High Loading Cases - P raax> > 

If the maximum loading applied to the plate is less than p^, then 
the actual variation of the loading can be analyzed directly. This is 
due to the fact that the hinge circle remains at the origin, and a single 
regime (AB) applies to the entire plate. For higher values of the 
loading, however, the hinge circle exists away from the origin and the 
actual shape of the loading is of prime importance. While the term 
"time varying loading" can refer, in general, to any variation with time, 
this term is restricted herein to loadings containing a single peak, that 
is, an initial increase in loading followed by a decay toward zero. This 
restriction is applied for several reasons. As mentioned previously, the 
main problems of interest to dynamic plasticity for circular plates are 
impact and ballistic problems, in which the loading is expected to have 
a single peak. In fact, estimation of the contact force experimentally 
(ref. ^1) by photographic means for a ballistic problem indicated that 
the force varied as 


F = F 0 t 2 e" pt (32) 

In addition, an excellent discussion of theoretical analyses of such 
ballistic problems (ref. 39) indicated that the method developed by 
Karas (ref. k2) also yields reasonable results for ballistic problems. 
Essentially, this approach assumed that a simply supported rectangular 
elastic plate was impacted by a spherical elastic projectile. Considering 
a Hertzian type of elastic impact, a difficult integral equation was 
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derived for the contact force history between the two elements and solved 
numerically. The numerical results for a specific case (20 cm by 20 cm 
steel plate with a thickness of 0.8 cm impacted by a 2 cm radius steel 
sphere) with a very low 1 meter/sec velocity, showed a response similar 
to equation (32). The magnitude of the maximum applied force was high 
(320 lb) and the duration of the loading was very small (12.8 4-sec). 
These values are quite surprising, since the velocity and mass of the 
projectile are small, and the resultant loading on the plate appears 
large. These results lend support to the often-applied assumptions 
that the loading starts at a maximum value and decreases rapidly. In 
any case, the theoretical results agree with the experimental results in 
predicting at most a single maximum peak. A short general discussion of 
both the increasing and decreasing phases of loading is given in the 
next two sections . 

Increasing loading above p^.- In this case, the plate has initial 
conditions determined from the solution with the hinge circle fixed at 
the origin as mentioned in the section entitled Low Loading Cases. 
Increasing the loading above p^ causes the hinge circle to move 
outward on the plate in a continuous manner. If the loading is allowed 

a discontinuity in time (such as an impulsive jump in loading), then all 

( <W 

except w and as required by continuity 

of the plate j can have discontinuities. The initial location of the 
hinge circle (p c ) can, therefore, be nonzero for a large initial impulse 
in loading. However, once the hinge has moved away from the origin, it 


may not have a discontinuous jump in position thereafter. This fact can 
be shown directly from the jump conditions derived previously. 

Therefore, it might be helpful to assume that as the loading 
increases, the hinge would move outward at a calculable rate, until it 
reaches a maximum value, after which it monotonically decreases toward 
the origin. Whether this is actually the case or not depends on the 
shape of both the radial and time variations of loading. There is the 
possibility that the shapes of these loadings may cause the hinge to 
move outward and inward several times before reaching the origin. If 
this is the case, great difficulty arises, for example, in determination 
of deflections, since the velocity expressions in the different regimes 
must be integrated for as long a time as they apply to any particular 
point. Finally, the yield surface restrictions on curvature rates 
(greater than or equal to zero in regime A) may preclude solution for 
some loading conditions. A discussion of this possibility follows 
equation (^7). 

Decreasing loading below •- In "this case, many of the 

accompanying difficulties of the previous section apply. As the loading 
is reduced, in general, the hinge may move inward or outward, and its 
exact location is important. For the case of a loading that is suddenly 
removed (see ref. 2), the hinge circle has been shown to decrease mono- 
tonically to the origin. If the loading is removed relatively 
quickly, then, it is reasonable to expect that the hinge will still move 
continually toward the origin. However, only careful solution of the 
hinge circle movement equation will tell how the hinge circle moves. 


BASIC ASSUMPTIONS 


In the solution of any general problem in plasticity a number of 
assumptions are needed to allow analytical solutions to be written. Some 
of the assumptions stated in this section have been discussed in earlier 
sections and will only be mentioned briefly for completeness. General 
discussion of these basic assumptions is available in references 1 
through 6. There are three general types of assumptions necessary for 
solution: assumptions on constitutive equations, assumptions on deforma- 

tions to be allowed, and geometrical assumptions required for the problem 
of interest. 

The material of the plate is assumed io be rigid-perfectly plastic. 
Rigid-perfectly plastic theory disallows the inclusion of work-hardening 
effects. Inclusion of work-hardening would necessitate a continually 
changing yield surface to be allowed, and the conditions for yielding 
would soon vary with the location on the plate and time, creating 
severe difficulties for solution. The material of the plate is also 
assumed to follow the Tresca yield criterion, which is a maximum shearing 
stress criterion (see fig. 4). Considering the boundary conditions to be 
only simply supported, the two regimes that apply are A and AB (see 
table 1), and the flow rules associated with these regimes define the 
strain-rate ratios from the normals to the yield surface as 

Regime A k T :k Q = 1 - p:p (0 < 0 < l) 

. . (33) 

Regime AB e r :ee = 0:1 

These flow rules allow solution for the deflections in many cases and are 
often applied most effectively by means of curvature rates. Assuming that 


55 


36 


normals to the original middle surface of the plate remain normal to the 
deformed middle surface, the strain rates throughout the plate depth are 
directly related to the curvature rates of the middle surface. In fact, 
the strain rates through the depth are 


. 3 u cPw , . u 

e= — = -z and e fi = - 

r Sr fcSfc 9 r 


1 b 2 v 


r drdt 


and the curvature rates of the middle surface are 


b 2 w , 1 dw 

— - and kq — — r-- 
ftr2 0 r dr 


so that 

e r :e Q = k r :k Q (36) 

It is also seen that the stress at each value of z must be at the yield 
condition for the plate to bend. Therefore, the stress above the middle 
surface of the plate is defined by one point on the yield surface, and 
the stress below the middle surface of the plate is defined by the 
diametrically opposite point on the yield surface. In this manner, the 
moments and stresses are related as 


n+ h 

r° r 

/ cr_z dz = a r 

/ z dz + a. 

J -h 

J-h 


This stress distribution is discontinuous at the middle surface, but such 
a discontinuity is allowable for rigid-plastic theory. 
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The plate is also assumed to be axisymmetric with respect to both 
loading and boundary conditions. This is a usual assumption for circular 
plate problems, and leads to the very desirable result of removing the 

shear stresses T r0 and t 0z identically. The twisting moments Mj. 0 

and the shear force Q 0 resulting from these stresses likewise vanish. 

If the plate is considered thin, then the vertical stress a z and the 
shearing stress r rz can be considered, on the whole, to be sm all in 

comparison to the bending stresses a r and <jq. This is the usual bend- 

ing theory for plastic plates and should be applicable for initial deflec- 
tions under dynamic loading. If membrane forces were allowed, then as 
soon as deflection started, additional load would be required to continue 
the deflection and, instead of allowing only bending moments to affect the 
yield surface, the inplane forces N r and Nq would have to be included. A 
four-dimensional yield surface would then be necessary with an appro- 
priate theory for interaction between these quantities at yield (see 
ref. 11). Such a generalization herein would obscure the desired insight 
into general plastic plate problems and is not included. Similarly, 
shearing deformation of the plate is neglected. This assumption is valid 
if the loaded surface of the plate extends over a region of the plate 
that is large compared to the thickness. In view of the allowance herein 
of rather general radial load distributions, the application to concen- 
trated loads is not made here. One final assumption is also made, that 
the loading on the surface of the plate is separable into radial and 
time functions, so that the influence of both effects can be Investigated 
separately. With these basic assumptions, then, the analysis can pro- 
ceed directly to evaluate the influence of load distribution. 


ANALYSIS 


Under the basic assumptions of the previous section, a typical plate 
element is shown in figure 6 with the applied forces and moments. In 
addition to the principal bending moments and shear force, the plate has 
a general surface loading q(r)p(t) and a resisting Inertia force as 
well. 


Basic Equations 

From the plate element shown In figure 6, the bt3ic differential 
equations can be derived. Taking the summation of forces and moments, 
and neglecting higher order differentials, these equations became 


^ (rQ r ) + rq(r)p(t) * pr ~ 


^ (m,,) - 


(38) 


where u represents the mass density per unit area of the plate middle 
surface. The shear force is eliminated from these two equations to yield 
one differential equation for the moments. 

~ (rMj.) - M q * -f |q(r)p(t) - n ^|jr dr (39) 

For a simply supported plate, the moments Mj. and Mg are equal at the 
center of the plate, and at the boundary the moment M r is zero. There- 
fore, from the Tresca yield condition (see fig. the center of the plate 
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corresponds to point A on the yield surface, and the edge of the plate 
corresponds to the point B. Since B cannot be a plasticity regime 
(table l), at most a central region of regime A and an outer annular 
region of regime AB exist, separated by a hinge circle located at 
r = p(t). The differential equation (39) must then be appropriately 
solved in each regime. The moment Mg is constant between points A and 
B of the Tresca yield hexagon, and consequently, it has been replaced by 
M 0 in the following sections. 


Solution - Regime A 

In this regime, M r = Mq and the differential equation (39) reduces 
to 



q(r)r dr - p 




dr 


0 


(40) 


To solve this equation for the "deflection rate" dw/dt, a separation of 
variables solution is written, so that 


dw 

dt 


q(r)R(r)T(t) 


(* 1 > 


Substituting this expression into equation (4o) and separating the radial 
and time variables yields 


ggili , C 
p(t) dt ' 


/; 


r dr 


J q(r)R(r) 


(* 2 ) 


r/r dr 


where C is a constant of separation. Solving these equations for the 
T(t) and R(r) expressions, 


... 


4o 


T (t) = £ J p (t)dt 


(J»3) 



dr 



m 


Taking the derivative of equation (44) with respect to r (or by inspec- 
tion), R(r) is seen to be simple i. Substitution of this value and 
equation (43), both into equation (4l), the deflection rate in this regime 
is found to be 


^ » a(£l r p(t )dt (45) 

bt 4 J 

Since the moments are already known, the solution is then complete. 

If the deflections are desired, then equation (45) must be integrated 
with respect to time. If the strain rates are desired, equation (45) 
must be differentiated with respect to r. Equation (45) shows the 
important fact that the velocity distribution in the interior region 
follows the radial distribution of the loading. This fact has generally 
been listed as an assumption previously, although only uniform radial load 
distributions have been investigated extensively. 

This expression has several restrictions on its applicability } for 
example, it applies only for r < p(t). If t^'e loading increases in time 
to a peak and then decreases, the equation (45) does not apply until 
P > Pj_ (defined in a previous section), since below Pj^ regime A does 
not exist. When p > p^, the regime A exists until the hinge circle 
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shrinks back to a point at the origin. This time must be determined 
directly from solution of the equation for hinge circle motion. Both 
before and after regime A appears, the hinge remains at the origin for 
a finite length of time. Otherwise, the plate remains rigid. Thus, the 
limits of integration in equation (45) are not directly defined. This 
is the basic difficulty caused by hinge circle movement. 

The flow rule also must apply in regime A. From equations (55) and 
(55), the flow rule in this regime states that 




dr2dt 




and 


K e “ - - 


a*w 


r brbt 


(46) 




Since p(t) and z r) are considered only positive quantities, then in 
regime A the following restrictions on radial load distribution must 
apply: 


d 2 q(r) 

dr^ 


< 0 


and 


dq(r) 

dr 


< 0 


(47) 


If the loading under investigation does not meet these restrictions, 
then two possibilities exist. Either regime A does not ever exist in 
that problem, or the present approach is inapplicable. It appears most 
logical to assume the plate remains entirely in regime AB, and for sane 
cases that have been analyzed, this was the case. Finally, it must be 
mentioned that the jump conditions between regimes A and AB across the 
hinge circle p must still be satisfied. 
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Solution - Regime AB 

Obtaining the solution in this regime is easier if the velocity- 
expression dw/dt is found first. Then, the acceleration tem can be 
evaluated in the moment equation, and the moment Mj. can be calculated 
directly. To determine the velocity distribution, the flow rule must be 
applied. In this regime, this is simply 


it r = - 


d^w 

dr^dt 


= 0 


(48) 


Solution of this equation for the velocity can be accomplished directly 
by separation of variables. With consideration for the requirement that 
the solution for velocity must be continuous across the hinge circle 
r = p(t), the general solution will be written instead as 


§=■ = q(p)ii(r,p,t) (4 9 ) 

Substitution of this equation into equation (48) yields a differential 
equation for f^(r,p,t) with the solution 

f^r^t) * rC x (p,t) + C 2 (p,t) (50) 

Realizing that the velocity must be zero at the supports (r * b), the 
functions C^(p,t) and Cg^jt) must be related as follows: 


C 2 (p,t) * -bC 1 (p,t) 


(51) 


^3 

Nov, the velocity must be continuous across the hinge circle r » p(t), 
and the velocities in both regimes are 
Regime A 




Regime AB 


( 52 ) 


bw 


l(p)(r - b) Cl (p,t) 


Equating these two velocity expressions at r = p yields C^(p,t) as 


Cl(o ’ t) s 57^1 / p(t)at 


' (53) 


Substitution of equations (50), (51) , and (53) back into equation (49) 
yields, finally, the velocity distribution in regime AB 


s - fer 1 / p(t)dt 


(54) 


This velocity expression for regime AB only applies for r greater than 
p(t). It is seen to be linear in r, so that the outer annulus of the 
plate deforms into a conical shape even in the general case* When the 
initial loading on the plate is less than p 0 , the plate remains rigid. 
Otherwise, regime AB applies for some region on the plate until motion 
ceases. The strain rate and deflection can be calculated directly from 
equation (54) by differentiation and integration, and the solution can 
be completed in this regime by determination of the moment Mj,. The 
expression for the acceleration tern d^/dt 2 is 


me et 
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and substitution of this expression into equation (39) yields the differ- 
ential equation to solve for M_. Note the rate of hinge circle motion 



*5 


l2(P> r ) 


Il(p» r ) =• f <l(05(r - 5)dt 
J P 

■f T (b - S)6(r - l)d| = ^(r 2 - p 2 ) 


(r + b) , , . (r 1 * - p^) 

(r5 - p5) + jj 


Two functions still remain undefined, Cj(p,t) and dp/dt. However, 
the moment M r has two boundary conditions to be satisfied, and these 


Mj * 0 and MJ = M q 
lr=b lr=p 


If it is realized from equations (60) and (6l) that when r =■ p both 
integrals are zero, the second of these boundary conditions yields 


C 3 (p,t) 


C 5 (p,t) =• pM Q 


Application of the first of boundary conditions (62) yields the differ- 
ent i:l equation to solve for p(t). This equation is discussed in the 
next section and will not be discussed here. However, the expression 
for Mj. contains dp/dt and, using the dp/dt equation, this term can 
be eliminated. After much algebraic manipulation, the moment expression 
can be written 


bM^^r) p( t) I I 2 (p,r) 

‘ V l(0,r) ‘ Il(0,b) 


pi8i! V- ^ ; ’ \ -VS? 
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Therefore, once the location of the hinge circle is known, the moment can 
be calculated, as well as^the velocity from equation (54), and the solu- 
tion for regime AB is complete. 

Hinge Circle Movement 

Consideration of the manner in which the hinge circle moves across 
the plate is a prime requirement for the solution of dynamic plasticity 
problems for thin plates. The rate at which the hinge circle moves aui 
be defined from the moment expression in regime AB, as mentioned 
previously. That is, equation ( 56 ) for the radial moment (with 
Cj(p,t) =■ pMq) jt be set equal to zero at the supports (r = b). This 
expression is taen solved for the hinge velocity dp/dt and yields 


where 


and 


dp bMp - p(t)pj 



= Ii(p>b) + 7g(p,t) 

j\(p) q(p)) I 2 (P,b) 

p 2 = <Tp-- (r^bT/r^by 


(65) 


Now, the basic manner of hinge circle movement can be discussed. 

The basic question of whether the hinge circle moves inward or out- 
ward depends on the sign of dp/dt. The time functions shown are all 
positive, but the p functions can have either sign and can cause the 
hinge to move in any direction. The integrals Ii(p,b) and Ig(p,b) 




( 66 ) 



The four quantities shown above are all positive, but their magnitudes 
may cause the hinge to move inward or outwa.-d depending on the shape of 
the q(r) function selected. The only possibility that exists in the 
general case is to use a computer to solve for p(t). A computer program 
to solve equation (65) for any general functions q(r) and p(t) has 
been written and a listing of the program is presented in the appendix. 

Exact solution of equation (65) is possible in only two general 
cases. The variables in the equation are directly separable whenever 
either p(t) or pj are zero or constant. When no loading is applied 
to the plate, p(t) is always zero and a trivial case exists, but when 
p(t) becomes zero after a fi*. .te length of time, then the motion of the 
hinge circle can be exactly calculated from that point onward. The 
function p(t) could also be zero in one other very important case. 

That case is any general Impulsive loading, in which p(t) is defined 
t>s zero, but the integral of p(t) is defined as a finite impulse applied 
to the plate. This general class of problems is very important , and a 
later section is devoted to the impulsive loading problem. The condition 
that p(t) is a constant occurs if the load is applied as p(> p 0 ) and 
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remains on the plate. Then, a fixed hinge circle problem is to be solved 
and equation (65) is not needed. 

The case in which p^_ is a constant can be determined directly as 
follows. The function q(p) must be found that makes the following 
equation true: 

Pi - / b q( 5 ) 6 (b - 6)dS - 4 r (b - p) 2 (b + 3 P) » C (68) 


Taking the derivative of this equation with respect to p and combining 
terms yields 


/. x dq(p) 

(b - p) — * -q 

dp 


Solution of this equation shows that the distribution under discussion is 

q(r) = C(b - r) ( 70 ) 

This is a triangular loading case and would allow the differential equa- 
tion on p to be separable. However, for this loading distribution, 
the hinge must only occur at the origin, and equation (65) is Inapplicable 
anyway. The deflection of the plate is conical in both regimes, and the 
definition of the hinf 1 circle breaks down in this case. Due to all these 
difficulties, the triangular loading case is not considered herein 

Matching Solutions at the Hinge Circle 
The solutions on both sides of the hinge circle must satisfy the 
Jump conditions derived in a previous section for hinge circles. These 
Jump conditions can be written out and proven to be satisfied in general. 


h'j 
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in regime AB then A-r.tward moving hinge 


* ( " t, / t t<r> ( ffi! ;, b> f*** ) at {*fr V**} « 


Po < r < p 


where t(p) is the time at which the hinge circle has moved to its 
current location, and t(r) is the time at which the hinge circle first 
reached the point r. These functions, as well as the functions of p 
that occur under the integral sign, must be determined from the solution 
of the hinge circle movement (eq. (65)) as functions of time to yield 
deflections. Now, if the Jump conditions on radial derivativer are to be 
written at the hinge circle, the deflections at points adjacent to the 
hinge must be writt u. If the hinge nas only moved outward, then the 
deflection on the inside of the hinge circle is expressed as the last of 
equations (71), and the deflection on the outside of the hinge circle is 
the second of equations (71) • If the hinge only moved inward, then the 
other two of equations (71) would apply, tor an outward moving hinge, 
the Jump in dw/dr when passing from the inside to the outside of the 
hlxge is 


1>J J t Q Wp - b) 



4 ( P )/p(t)dt .niriajd r* t)at 



dt 


( 72 ) 
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If r = p all the integral terns cancel directly, and all dt /dr terms 
(when multiplied by dp/dt) also cancel, satisfying the Jump condition 
(l4). For the other case of an inward moving hio^e, the same equa- 
tion (72) is identically found. If the hinge moved inward and outward 
several times, then each time it passed a given point r another inte- 
gral factor would have to be added to the deflection expressions to 
denote the change in regime. However, the Jump condition on Sw/Sr is 
directly satisfied whether the hinge moves either inward or outward.. 
Also, the Jump condition on S'Tr/Sr 2 can be proven, in general, by 
carefully differentiating equation (72) and letting r = p for both 
inward and outward moving hinges. 

The Jump condition of equation (15) can be verified, in the general 
case of a moving hinge circle, by means of the velocity expressions 
already derived in both regimes. Consider the Jump as the magnitude 
when passing from the inside to the outside of the hinge circle, then 


and 


d 2 * _ Sq(r) f p(t)dt _ q(p) f p(t)dt 

SrStJ Sr J 4 (p - b ) J ^ 


lr=p 


(75) 


■ / p(t)dt {(77 


-q(p) 


1 dq(p)) dp 
b)^ (p - b) So j dt 


( 74 ) 


r»p 
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How, multiply equation (73) by dp/dt and add to equation (74), letting 
r * p, so that the jump condition becomes 


q(p) 


p(t)- 


q(p)(p - b) 
p(p -b) 


p(t) - 


(p-b) 




p(t)dt q(p) 





By close inspection of equation (75), it can be seen that all terms are 
canceled by similar terms and this jimp condition is identically 
satisfied. 

The jimp condition cm moment Kp can also be easily 'verified in 
the ge ner al case. Since Up = Mq in regime A, the jump condition 
(eq. (16) ) defines the value of the derivatives of Mp in regime AB as 
the hinge circle is approached. Then, in regime AB (eq. (64)) 


Wo /l dlgtar) l2(p,r)\ ) i dl-^r) 

dr I 2 (p,b)\r dr " r2 / “ ~ 


and 

bt 


Il(p,r) 

r2 


T l(p> b )/ l dlg(p,r) 
Ig(p,b)\ r dr 



(76) 


f / X dlp(p,r) dlp(p,b] 

^Ig (p,b) gj ■ - Ig(P,r) 

(ia(p.b)> a 


P(t) dp 

~dt 


|dl 1 (p,r) 
dp 


d I 1 (p>b)I g (p,r)) idn(t)/ T . . Ij(p,b )^r)\ 

“dp ^(p,b) Jr dt |l (p,r) " Ig(p,b) J 


(77) 


§ 
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Noting that I 1 (p,r) = I 2 (p,r) * 0 as r approaches p, and that adding 
equation (77) to dp/dt times equation ( 76 ) should yield zero (the jump 
condition), then 


Wq ^Igtpjr) 
pl 2 (p,b) 5r 


r*p 


dp _ p(t) dli(p,r) 
dt p 5r 


dp 

dt 


r=p 


+ p(t) 


I^(p,b) dl 2 (p,r) 


p!g(P»b) dr 


r=p 


do »o **g(p.r> 

dt pl^pjb) dr 


lr=p 


dp 

dt 


p( t )/d I i(p, r ) 


dp 


I^Pib) dlg(p,r) 


r=p 


Ig(p»b) dp 


f£-0 

dt 


(78) 


The derivatives that new renal” can be evaluated directly, and 


dl.(p,r) nr 

-V~ -J f * (t >‘ « 

» | (r2 - p2) - i (r3 - o5) 


dl x (p,r) 

dp 


dl 2 (p,r) 

§p 


-q(p)p(r - p) 


(79) 


= -p(r - p)(b - p) 


Since one of these Integrals appears In each term of equation ( 78 ) and 
they all became zero as r approaches p, the Jump condition on moment 
Mj. is identically satisfied. 

Tnitiai Hin ge Circle Location - Impulsive Loading 
Definition of the motion of the hinge circle is the primary diffi- 
culty in dynamic plasticity problems. Another basic difficulty is the 
definition of the Initial location of this hinge circle. If the loading 
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varies from a low value up to a high peak value and then dawn again, the 
hinge obviously forms at the origin as p = p Q , and moves away from the 
origin at p = p^. It also reaches a maximum value at some later time 
and decreases to zero thereafter. This type of smooth motion can be 
analyzed directly without difficulty. However, if a large step discon- 
tinuity in time occurs initially in the loading (as in Impulsive loading 
and other problems), the hinge may form at a location away from the origin. 
This is the only case in which the hinge may have such a "jump, " but this 
case is important, and location of the initial hinge circle is required. 
For regime A to exist in the center of the plate, the curvature zates ic r 
and kq must be non-negative. Therefore, the initial hinge location p Q 
must be the smaller of the two values that make K r and kq non- 
negative in regime A. Therefore, 




a 2 q(r) 

dr2 


> 0 


1 d^w 


* 0 " ‘ r drdt 



> 0 


(80) 


The only restrictions, then, on initial hinge circle location are that 
p 0 must be the smaller of the two values 


dq(p Q ) 

8p 0 


< 0 


and 


^q(p 0 ) 


(81) 


The shape of the loadings is seen to be very important for large initial 
step loadings. 
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If the loading is, indeed, a totally impulsive loading, then one 
additional restriction on p Q is possible. For impulsive loading, the 
load has been applied to the plate instantaneously, and the energy 
imparted to the plate must be dissipated in plastic deformation as the 
plate stops. If the hinge circle stayed at p Q or increased in size, 
the energy of plastic deformation would not increase as it should. There- 
fore, for impulsive loading the hinge circle must decrease or dp/dt must 
be less than zero. For an impulse, p(t) = 0, and J p(t)dt = (iV 0 is 
the impulse applied to the plate. In this case, then, the initial hinge 
velocity becomes 


dp 

dt 


bMp 

P=Po P v 0 p 2 


-bMo 12 


uv 0 (b + 5p 0 )(p 0 - b)' 


2 pq(po) . <l(Po) 


t 


(b - P 0 > 




(82) 


This velocity is less than zero whenever the following inequality is also 
true: 


dq(p Q ) q(p Q ) 

ap o (b - P 0 ) 


( 85 ) 


Thus, in summary, if a large Initial discontinuity in load occurs, 
the smallest value of p Q determined from equations (8l) applies. If 
the loading is truly impulsive at t = 0, then the lowest p Q from 
equations (8l) or equation (83) must be used. The magnitude of the load- 
ing has no effect on these expressions, only the shape of the curve q(r). 


56 


Exact Solution - General Impulsive Loadings 
For a general impulsive loading, the equation for the motion of the 
hinge circle can be solved exactly, and since most previous work har been 
concerned with this type of loading, this general solution is presented 
here. For the impulsive case p(t) = C and J p(t)dt = P v o so that the 
equation ( 65 ) becomes 


bMo 

Po dp * dt (84) 

P V o 

and, since the variables are separated, direct integration yields the 
solution. Using integration by parts, the p 2 integral becomes 

/**-/ I * (p ' b) ^((7?ry) d0 - (7TT) 

+ J pq(p)(p - b)dp ( 85 ) 


and after integration of equation (84) 




q(p)(b - p) 2 (b + 3 p) 


12 


f q(p)(p 2 - bp) dp (86) 

J Po 


Po 


Consideration of the initial conditions (z = 0, p = p Q ) yields =• 0. 
This is the general solution for hinge circle motion. Note that the only 
integrations that remain are J' p^q(p)dp and J pq(p)dp, and these 
should be directly obtainable for nearly any function. The initial hinge 
location p Q is found to be the smallest value for either equations (8l) 
or equation (83). 
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Also, the velocity distributions become (from equations (45) and 


(5*0) 


Regime A 
Regime AB 




8w 

5t 


» q(p) 


(r - b) 


(87) 


The moment Mj. in regime A is M Q , and in regime AB the moment Mj. 
becomes directly from equation (64) 


Mr » Mo 


12bMo 


r(b - p) 5 (b + 3p) 


. p2) . <2LL±>< r 3 . p3) + 


( 88 ) 


The hinge circle reaches the origin at the time t^, defined by the left- 
hand side of equation (86) with p = 0 on the right-hand side. 


tl 


4V 0 b c 

12Mq 


jjq(o) + q(p 0 )(l - t) (*■ + 5 Tr) + TT 


q(p) 



After the hinge reaches the origin, — =■ 0 and p = 0 and since 

at 

p(t) = 0, the acceleration in regime AB vanishes (see eq. (55) )> and the 
velocity distribution in regime AB (eq. ( 87 )) becomes 


8w 

dt 




(90) 


For values of t > t^, this velocity cannot remain constant in time 
since energy is still being dissipated. Following the same approach as 
outlined in reference 5 for a specific q(r), the velocity distribution 
is assumed to decrease linearly with time (t > t^) until the plate stops. 
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The velocity distribution then becomes 


8w 8w 
8t 8t 


+ ifr* (ti - t)(l - £) 
b H \ W 


(91) 


(t > t x ) 

The plate finally comes to rest at the time t* when the velocity drops 
to zero, so that 


„V 0 ,(o)b‘ 
1 12M n 


(92) 


The solution for M r in regime AB (eq. (88)) for p = 0 reduces to 


^-- 4 - 44 ) 


(93) 


Thus, it is noted in this case that the shape of the loading affects only 
the solution up until the time t^, and the final conical deflection that 
decays to zero velocity is the same regardless of the loading shape. 


Determination of Final Results 

After selection of the loading functions q(r) and p(t), the 
solution of a general plastic plate problem still remains a complex 
operation. First of all, the loading must be checked to be certain 
p > p D and a hinge will form. Also, p > p-^ should be checked to 
determine if the hinge occurs away from the origin. Assuming it does, 
the dp/dt equation must be solved. TLc program in the appendix can be 
used for general cases simply by adding the required definition of the 
functions in the proper subroutines. After the variation of p with time 
is determined, this value of p must be substituted in the proper 
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expressions for velocity and moment to find their distributions at that 
particular time. If other quantities are of interest (strains and strain 
rates, for example), then the velocity expressions (eqs. (45) and (54)) 
and deflection egressions (such as eqs. (71)) must be integrated or 
differentiated. In general pro 1 :* sms, these manipulations become very 
difficult if the motion of the hinge circle is more complex than simply 
an outward then an inward movement. In the simple case of impulsive 
loading, the exact expressions can be written in a straightforward manner. 

Example Cases 

To illustrate the influence of the time variation of loading, several 
example cases were computed for general loading functions. To allow a 
general type of loading that may be expected in practice, the loading 
q(r) was selected to be Gaussian, and the time variation of the load was 
considered to be essentially exponentially decaying. Since the combina- 
tion of general radial and time functions has not been attempted before, 
some latitude exists as to the exact shapes to select. For the examples 
herein, then, the radial distribution is 

q(r) = e" a2r2 = e _r2 (9*0 

This load decreases from a maximum at the origin to approximately 58 per- 
cent of it 8 maximum value at the support (r = b = l), and is non-uniform. 
The time function was selected so that a uniform load was applied until 
a small total impulse had been given to the plate, then an exponential 
decay was assumed. The general shape considered is shown in figure 7* 

The loading remained constant until an Impulse of 1 x 10~^ had been 
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applied. Thereafter, the loading dropped of exponentially as e* 0 *, and 
the value of a decided hov rapidly the loading decreased. The total 
impulse applied to the plate during the unloading phase was also held 
constant at one. Attention is centered on solution for the hinge circle 
movement, since this motion is considered the basic kernel of the problem. 

The first sample case is shown for a = 1 in figure 8. Since these 
cases all have an initial jump in loading, the hinge starts at a location 
other than the origin. In this case, the hinge starts at p Q = /I/2 
(as required by the radial load distribution), but it quickly moves to 
the origin of the plate. This same type of motion has been seen pre- 
viously for uniform loading and is an example of a continually decreasing 
hinge circle. 

The second sample case is shown for a =* 1 x 1(P in figure 9. In 
this case, the first example of an outward moving hinge circle is seen. 

Due to the radial load distribution again, this hinge also starts at 
p Q = J2/2, but soon after the load variation starts, the hinge moves out 
to a maximum value of 0.98^ before it slowly moves inward. Note that the 
hinge circle moved outward even though the loading is decreasing with time. 
Previous results on very high impulsive loadings (ref. 2) indicate that 
the Initial hinge circle location was seen to approach b r as the 
loading approached infinity. Therefore, this example of a high initial 
loading seems to verify this effect. Note also that even after the same 
total impulse has been applied in both cases, the location of the hinge 
circle is vastly different from case 1 due to the time factor only. 

Perzyna (ref. 37) , for a uniform radial distribution and different time 
expressions, states that the influence of time variation of the loading 
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is unimportant in determination of hinge circle location, but this is 
not generally the case. The velocity and moment distributions have been 
calculated for case 1 and are shown in figures 10 and 11, respectively. 

As expected, figure 10 shows the velocity distribution is linear (or 
conical.) in the region outside of the hinge circle. The shape of the 
velocity distribution inside of the hinge circle is non-unifoim and of 
the same general shape as the loading function q(r) . However, it can 
be seen that the velocity in the central region of the plate is contin- 
ually increasing. This effect was not seen previously (ref. 3 6 ) for 
impulsive loading of this general shape, and is caused by the time varia- 
tion of the loading. The moment distributions in figure 11 are seen to 
smoothly decay from M 0 at the hinge circle to zero at the support. 

There appears to be no discontinuity in slope at the hinge, and the 
curves tend to have a slight reverse curvature near the support (especi- 
ally for later times). 

The velocity and moment distributions lave also been calculated for 
case 2, and are shown in figures 12 and 13, respectively. Since the 
hinge circle has been allowed to move outward, the velocity distributions 
in the center region of the plate have a reverse curvature in them at 
r » 0.707 b. This same reverse curvature exists ir q(r), and the shapes 
of the velocity distributions agree with q(r) out to the hinge circle. 
These distributions are also conical in the outer region of the plate. 

The increase of velocity with time in the central region is even more 
pronounced in this case than in case 1. The moment distributions in the 
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plate for case 2 (fig. 13) show that for early tines the moment distribu- 
tion is similar in form to case 1. However, as the hinge moves outward 
with time, the moment distribution becomes steeper and the reverse curva- 
ture seen only slightly in figure 11 becomes very pronounced in this 


case. 


CONCmSIONS 


From the results and equations shown herein, several important con- 
clusions are evident. The equations derived here considering bending 
deformations only are seen to be more general in form than existing solu- 
tions, and reduction to the existing cases is direct. For example if the 
loading is considered uniform in r and impulsive or step-wise uniform 
in time, the equations derived directly for such cases by Hopkins and 
Prager and Wang (refs. 2 and 5) appear exactly. Also, if the radial 
load distribution is considered uniform, and a general function of time 
is allowed (but assuming only inward hinge circle movement), the non- 
linear equations of Perzyna (ref. 37) are found exactly. The conclusion 
of Perzyna that time variation is unimportant appears to be caused by 
an unfortunate choice of example time functions. He solves the specific 
non-linear equations for his example, and does not present any means 
for evaluation of his numerical method of solution. 

If the loading on the plate is considered to be a distributed 
Gaussian loading in r and impulsively applied, the equations derived 
directly for this case by Thomson (ref. 36 ) appear exactly herein. 

These two papers (by Perzyna and Thomson) are the only two papers avail- 
able at present that allow variations of the loading, one in r and 
the other in t, and both sets of equations are included in the general 
expressions herein. In fact, the solutions currently available for 
bending theory are found to exist as special cases of these general 
equations. 
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Two other results presented herein are also of direct interest. The 
equations for impulsive loading of general radial shape are shown to be 
directly solvable, and expressions for the solution are written. Many 
cases of practical interest are, therefore, seen to be available now in 
a direct manner without numerical integration or other approximation. 

Also, the general load-carrying capacity of circular plates under arbitrary 
radial load distribution is presented for the Tresca yield criterion. 

This load value (p Q ) is known to yield a lower (conservative) bound for 
the load-carrying capacity and can be used effectively for design purposes. 
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APPENDIX 


GENERAL COMPUTER PROGRAM FOR HINGE CIRCLE MOTION 

A general computer program that solves the dp/dt equation 
(eq. (65)) for arbitrary functions of r and t has been written in 
Fortran IV computer language, said a listing of the deck is included in 
this appendix. Many comment cards have been included in the deck to 
directly define the expressions and quantities involved and should be 
essentially self-explanatory. If a new radial distribution is to be 
investigated, it must be defined in the function Q(R), and its first 
and second derivatives written as the functions QP(R) and QFP(r), 
respectively. If the exact expression for I^(p,b) can be written, 
this expression should be put into the subroutine I 1 EX. If it cannot 
be expressly written, the main program will numerically integrate 
(simply set IEX = 0 ). Similarly, if a new time function is to be 
investigated, it must be defined in the function P(T). If the exact 
expression for J P(t)dt can be written, this expression should also 
be included in the subroutine FTEX. If not, the main program will 
integrate numerically (simply set NEX = 0 ). A sample case is also 
included, and the initial printout shows the procedure used in the 
machine to determine the initial hinge circle location. If a totally 
impulsive loading is applied to the plate, an additional possible 
restriction on p 0 (eq. (85)) must also be considered. However, for 
impulsive loading it is more reasonable to use the exact expressions 
derivable as shown in the text. 
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PROGRAM GtNPLAS (INPUT, OUTPUT ) 

THIS IS THE MAIN PROGRAM FOR SOLUTION OF THE ORHU/DT EQUATION 
AND IS THE SAMt FOR ALL GENERAL TYPES OF LOADINGS 
S^RADIUS OF ThE PLATE ,MZ=INITIAL YIELD MOMENT M ZERO 
II, 12, RHOl, RHOZ ARE DEFINED FUNCTIONS OF THE HINGE LOCATION 
TZMNITIAL TIME, A SHORT TIME DURING WHICH THE LOAD IS CONSTANT 
OELT IS AN EXPECTED TIME STEP USED TO PLOT RHO VS TIME 
E 1 , E 2 , ARE RELATIVE ERRORS USED IN ITR? TO FIND RHO INITIAL 
M AX I = MAXIMUM NUMBER OF STEPS ALLOWED IN RHO VS TIME RESULTS 
RB.RE.Rt ARE THE STARTING .ENDING .AND INCREMENT VALUES OF RHOZ 

USED IN THE ITR? ROUTINE 

K1,K2,K3,K4 AkE INTERMEDIATE ANS. USED IN THE RU''G-<UTTA SOL. 

A IS A CONSTANT THAT TYPIFIES THE RADIAL LOAD VARIATION 

ALPHA IS A CONSTANT THAT TVPIFIES THE TIME LOAD VARIATION 

I EX=C USE NUMERICAL INTEGR. FOR I1,IEX=NON-0 EXACT EXPR. IS USED 

NEX=0 USE NUMERICAL INTEGR. FOR PT,NEX=NON-0 EXAtT EXPR. IS USED 

OERSUB IS THE ROUTINE THAT CALCULATES DERIVATIVE DRHO/DT 

P. HOI NT IS THE ROUTINE THAT CALCULATES THE INITIAL HINGE LOCATION 

ANS ALWAYS REFERS TU THE DERIVATIVE DRHO/OT 

X AND Y REFER TO TIME AND RHO IN RUNG-KUTTA SOLUTION 

ITR2 IS A MACHINE ROUTINE THAT SOLVES FOR THE ZEROS OF A FUNCTION 

mgauss is a machine routine that integrate? numerically 

DIMENSION RESULT l 2 J 
REAL K1,K2,K3,K4,MZ,I1,I2 
COMMON/BLKl/B.A, ALPHA 
COMMON /BLK2/MZ 
C0MM0N/BLK3/I 1 , I 2 , RHOl ,RH02 
COMMON/BLKA/tl ,E2,MAXI,TZ 
COMMON/BLK 5/IEX.NEX 
COMMON/ BLK6/ IRR 

NAMELIST/INPUT/TZ,MZ,B, A, DELT, ALPHA, MAXI ,E1,E 2,1 EX.NEX 

1 CALL DAYTIM (RESULT! 

PRINT 2, RESULT 

2 FORMAT(*1DATE*3XA10»5X*TIME*3XA10/* ROB I NSON-WE I DMAN , SRD-A2058, R 
10P-308*/* GENERAL PRESSURE LOAOING ON PLASTIC THIN PLATFS*//) 

READ INPUT 

IF ( IEX.FU.O) PRINT 25 

IF (NEX.Ev.O) PRINT 30 

25 FORMAT! *ONUMER IC AL INTEGR. USEO FOR II*) 

30 FORMAT ( *ONUMER I CAL INTEGR. USED FOR PT*) 

3 CALL RHOINT (RHOZ) 

PRINT 20, MZ, ALPHA, B, A, DELT, RHOZ, TZ 
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20 FORMAT ( 5X*MZ=*fc 1 5 . 8 , 5X* ALPHA**E1 5 . 8/ 

16X*B=*El5.6t9X*A**El5.8/ 

23X*DELT=*E15.b,6X*RH0Z=*El5.S/ 

35X*TZ =*E1 5 .8// ) 

PRINT A 

A FORMAT! 10X*T*» lt»X*RHO*» 16X*1 l*t 16X* 1 2* » 1 5 X*RH01 *» l AXARH02* * 1AXMNS 
l*//) 

R HO=RHOZ 
T = TZ 
H- DEL T 
I =0 
5 X-T 
Y = RHO 

CALL OERSUe ( X » Y i ANS) 

PRINT 8, X t Y , 1 1 , 1 2 > RHO 1 > RH02 » ANS 

8 F0RMAT(7!3XfclB.8)) 

IF (RHO.LE.C.O) Gl) TO 1 

K I =ANS*H 

X=T+H/2. 

Y-RH0«-Kl/2. 

CALL OERSUf* ( X » Y» ANS) 

K2=ANS*H 

Y=RH0+K2/2. 

CALL DERSUb IX, Y, ANS) 

K3 =ANS*H 
X = T+H 
Y---PH0+K3 

CALL OERSUb ( X » Y , ANS) 

KA=ANS*H 

DELRHO=!Kl*2.*K2+2.*K3+KA)/6. 

1 = 1 + 1 

IF ! I.GT.MAXi) GO TO 6 
IF(ABS(DELRH0).G£..15*B) GO TO 9 

IFIABSIDELRHO) . Lfc . . 00 1 *B. AND. ABS I A*IS ) . GE . 1 . E+ 2 ) GO TO 10 
T = T+H 

RHO=RHO+DfcLF.HO 
GO TO 5 
fc PRINT 7 

7 FORMAT! BX*Ma;UMUM NUMBER OF STEPS TAKEN*//) 

GO TO 1 

9 DEL T = DELT * . 2 
12 PRINT 1 5 * DcL T 

IS FORMAT! 5X*DtlT HAS SEEN CHANGED TO *E15.8) 
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H*DELT 
GO TO 5 

10 DELT«DELT*5. 
H-OELT 

PRINT 15 t DkLT 
GO TO 5 
STOP 
END 


nnoonn 
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sub; .-tine DtRsuB «t,rho*ans> 

THIS ROUT 1 Pit CALCULATES THE DERIVATIVE DRHO/DT 

ANO IS THE SaHE FOR ALL GENERAL TYPES OF LOADINGS 

FUNC1 * FUNL2 tFOFXI t FOFX2t SUHl tSUM2 ARE ALL ADOED FOR INTEGR * ROUT* 

IlEX CONTAINS THE EXACT EXPR* FOR II IF ANALYTICALLY KNOWN 

PTEX CONTAINS THE EXACT EXPR. FOR PT IF ANALYTICALLY KNOWN 

FOR OTHER DEFINITIONS OF FUNCTIONS SEE MAIN PROGRAM 

EXTERNAL FUNCI.FUNC2 

DIMENSION SUMMn,FOFXim*SUM2m.FOFX2Ill 

REAL I2vIifM£« NUM » 1 2 A 

COMMON/BLKl/B, A, ALPHA 

COMMON/ BLK2/MZ 

COMMON/BLK3/ 1 1 » 12 » RHQ1 9 RH02 

CQMMON/BLKS/ I EX, NEX 

I 2=-I RHO-tt ) **3*< 3.*RH0+B1/I2. 

!2A=I9-RHO*v|3.*RHQ+B>/12. 

IF < IEX.EU.O) GO TO l 
CALL I1FX Ul*RHO) 

GO TO 2 

1 INT=f B-RHGJCIO. 

CALL HGAUSS IkHQ, B, INT f SUNI , FUNC1*F0FX1 t 1 > 

Il = SUMim 

2 RH02= I 2A* { gP ( kHliJ * ( RHO-B )-QI RHO I ) 

RH01=Il+Q<RH0»*I2A*IRHQ-B) 

IF INEX.EU.O) GC TO 3 
CALL PTEX (PT t T) 

GO TO A 

3 CALL MGAUSS « 0 . ,T , 30,SUM2 ,FUNC2*FOFX2* I ) 

PT=SUM2U) 

A DEN=PT*RH02 

NUM=B*M2-M T ) *P.HU1 

ANS=NUM/DEN 

RETURN 

ENO 
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SUBROUTINE RHOINT tRHOZ) 

THIS ROUTINE CALCULATES THE INITIAL HINGE LOCATION 

AND IS THE SAME FOR ALL GENERAL TYPES OF LOADINGS 

ICODE IS AN ERROR DIAGNOSTIC WRITTEN TO DEFINE WHERE THE MIN. 

VALUE FOUND BY I TR2 IS LOCATED. 

I CODE=0 PROPER RETURN... ICODE=3 THEN RHO OUTSIDE OF THE PLATE 
ROIN IS ANOTHER EXPR. WRITTEN FOR RHO INITIAL. 

FOR OTHER DEFINITIONS OF FUNCTIONS SEE MAIN PROGRAM 

CCMMON/BLKI/B* A# ALPHA 

CCMMON/BLK<,/El,E2,HAXI ,TZ 

C0MM0N/BLK6/IRR 

EXTERNAL UPP 

EXTERNAL UP 

R B-. 001*B 

RE=l.O*B 

R I =. 1 *B 

CALL I TR2 IRHOZ, Rb ,RE , R I ,QPP , El, E2, MAXI, ICODEI 
IF (ICODt.NE.O) PRINT 1, ICODE 
PRINT 2 

2 FORMAT ( *0*// ) 

1 RR=0 

l FORMAT! *OEPRQR FROM ITR2- QPP IC0DE=*13/) 

IF! ICODE. EO. 3) RH0Z=0. 

I F I I COOF. EU.3. AND.UPP (O.).LE.O.) RHOZ=B 

PRINT 2 

RC=.001*B 

CALL ITR2 (ROIN,RC,KE, R I ,QP , E 1 ,E2 , MAX I , ICODE ) 

IF I I CODE . NE . 0 ) PRINT 3, ICODE 
PRINT 2 

3 FORMAT! *OERROR FROM ITR2- QP IC0DE=*I3/) 

IT! ICODE. EQ. 3) RGIN=0. 

IF ( I CODE. £U. 3. AND. UP (O.).LE.O.) ROIN=B 
PRINT 2 
I RR = 1 

IF(ROIN.LE.RHuZ) RHOZ=ROIN 

RETURN 

END 
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SUBROUTINE FUNGI <X,FOFXl) 

THIS ROUTINE IS SIMPLY AN INTEGR • NECESSITY 

ANO IS THE SAME FOR ALL GENERAL TYPES OF LOADINGS 

DIMENSION FOFXim 

C PMMON/6LK I /tt f A »AL PH A 

FOFX1C 1 >=<B-X)*Q(X)*X 

RETURN 

END 
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SUBROUTINt FUNC2 IX.FOFX2) 

C THIS ROUTINE IS SIMPLY AN INTEGR. NECESSITY 

C AND IS THE SAME FOR ALL GENERAL TYPES OF LUAOINGS 

DIMENSION FuFX2 ( l ) 

FOFX2 ( 1 ) =P< X) 

RETURN 

END 
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FUNCTION UtRi 

THIS ROUTINE IS WHERE THE RADIAL DISTRIBUTION IS DEFINED 
AND MUST bE INDIVIDUALLY WRITTEN FOR EACH LOADING SHAPE... 
C OMMON/BLKI/ B» A » ALPHA 
Q=EXP(-A**2*R*R) 

RETURN 

END 
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FUNCTION OPiR) 

C THIS ROUTINE CALCULATES THE FIRST DERIVATIVE OF THE RADIAL LOAD 

C DISTRIBUTION AND MUST BE WRITTEN FOR EACH LOADING SHAPE... 

C OMMON/8LKI/B , A * ALPHA 
C0MH0N/BLK6/IRR 
QP=-2.*A*A*R*tXP(-A**2*R**2l 
IF ( IRR.fcQ.O) PRINT l,R,QP 
l FORMAT! 5X*RHG=*fcl5.8, 5X*QP=*E1 5. 8 ) 

RETURN 

END 
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FUNCTION QPP(RHO) 

THIS ROUTINE CALCULATES THE SECONO DERIVATIVE OF THE RADIAL LOAD 
DISTRIBUTION AND MUST BE WRITTEN FOR EACH LOADING SHAPE... 
COMMON/BLKI/B, A, ALPHA 

QPP=EXP(-A**2*RH0**<:)*(-2.*A**2+4.*A**A*RH0**2> 

PRINT I, RHO.UPP 

1 FORMAT I 5X*RN0=*E 15 .8* 5X*QPP=*E 15. 8 1 
PETURN 
ENO 
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FUNCTION FIT) 

C THIS ROUTINE IS WHERE THE TIME DISTRIBUTION 

C AND MUST bE INDIVIDUALLY WRITTEN FOR EACH TIM 

COMMON/BLKl/b,A,ALPHA 
CONST =.lF-3/ ALPHA 
IF IT. LE. CONST) GO TO l 
P=ALPHA*EXPl-ALPHA*lT-CONST) ) 

RETURN 
1 P=ALPHA 
RETURN 
FND 


S DFFINEO 
VARIATION.... 
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SUBROUTINE 11 tX (Il.RHOI 

THIS ROUTINE CALCULATES EXACT EXPR. FOR II IF ANALYTICALLY KNOWN 
AND MAY BE WRITTEN UR OMITTFD FOR EACH SPECIFIC RADIAL LOADING 
TYPIFIED BY A PARAMETER A 

SET I EX=NUN-0 ON JINPUT CARD TO USE THIS ROUTINE 
REAL II 

COMMON/BLKI/bt A, ALPHA 
P I =T. lA159«!ob35C979 

I 1 = ( EXP ( - A**2 *KHO*''2 )* ( A*B-A*RHO) / 2. +SQR T( P I ) *( ERF I A*RHO ) - 
IEKF(A*B> >/A. )/A**3 
RETURN 
END 
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SUBROUTINE: PTEX (PT,T) 

THIS ROUTINE CALCULATES EXACT EXPR. FOR PT IF ANALYTICALLY KNOWN 
AND MAY BE WRITTEN UR OMITTFO FOR EACH SPECIFIC RADIAL LOADINC, 
TYPIFIED BY A PARAMETER ALPHA 

SET NEX=NUN-0 ON JINPUT CARO TO USE THIS ROUTINE 

COMMCN/BLK1/u,a, ALPHA 

X=ALPHA*T 

IF IX.GT.I.E-4J GO TO I 

R T = X 

RETURN 

1 PTrl.o*I.c-4-tXP(-(X-l.E-A) ) 

RETURN 

END 
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DATE 04/23/68 TIKE 13.26,25. 

ROBI NSON-WE I DM AN f SRD-A2Q58, RDP-308 

GENERAL PRESSURE LOADING ON PLASTIC THIN PLATES 


RHQ= i.COOOOOOOE-03 
RHG = 1 • 0 1QC00006-01 
RHO= 2.010000006-01 
RHO = 3. 010000006-01 
RHO = 4. OlOOOCOOE-Ol 
RHO= 5.01 OOOOOOE-Ol 
RHQ = 6.01000000E ’01 
RHO= 7.01000000E-01 
RHO = 8 • OlOOOOOOE-Ol 
RHO* 7.510000006-01 
RHO = 7.260000006-01 
RHO- 7. 1350000QE-01 
RHO= 7. 0725000 OE -01 
RHO = 7.041 25000E-01 
RHO* 7 • 056 875006-01 
RHO= 7 . 06468 750E-0 1 
7H0 = 7. 06 859375 6- 01 
RHO* 7.070546876-01 
RHO= 7.071 52 344E-01 
RHO= 7.07103 5 .66-01 
RHO= 7.071279306-01 
RHO= 7.0711 5 723E-01 
RHO* 7.071096196-01 


RHO* 0. 


RHO* 1 • 00000 000 E- 03 
RHO* 1 • 010000006-01 
RHO- 2.Q1000000E-01 
RHO* 3 OlOOOOOOE-Ol 
RHO* 4.010000006-01 
RHu= 5.01000000E-01 
RHO* 6.010000006-01 
RHO= 7.0 10000006-01 
RHO= 8.010000006-01 
RHO= 9.01000000E-01 
RHO* 1.000000006+00 


WPP--1. 999994006+00 
UPP=- 1.9393 11 83E+00 
QPP=-1. 765603316+00 
QPP*-1.49575084E+00 
QPP*-1. 155259556+00 
QPP =-7.749063956- 01 
QPP--3.86882181E-01 
QPP*-2. 104239646-02 
QPP* 2 . 98 1 82765E-01 
QPP* 1.4 564 7901E-0 1 
QPP * 6.393494546-02 
QPP = 2.183543426-02 
QPP* 4.9 1 34 1929E-04 
QPP=-1.02521363E-02 
QPP*-4. 874509846-03 
QPP=-2. 190107186-03 
QPP=-8. 490128136-04 
QPP --1. 787429156-04 
QPP = 1.563226496-04 
QPP— 1. 120434906-05 
QPP- 7. 256059616-05 
QPP = 3.067848516-05 
QPP = 9.737158416-06 


QPP=-2 • 0 00000 OOE +00 


QP* - 1 • 99999800E — 0 3 
QP — 1. 999498726-01 
UP — 3.86082504E-01 
QP=-5. 498560116-01 
UP — 6.828721196-01 
QP=-7.79577637E-01 
QP =-8 .376003826-0 l 
wP — 8.57699725E-01 
QP=-d. 433711456-01 
UP *-8*001917916 -01 
yP=-7. 357588826-01 


ERROR FROM ITR2- QP ICOUE* 
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TABLE 1— POSSIBLE PIASTICITY REGIMES FOR CIRCULAR PLATES 
WITH VARIOUS YIELD CONDITIONS 


Yield 

conditions 

Possible regimes 

Simply supported 

Clamped 

von Mises 

AC 

AE 

Tresca 

A, AB 

A, AB, HD, D 

Johansen 

A, AB 

A, AC, C, CD 

Haythornthwaite 

AB, BC 

AB, BD, DE 
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Figure 2 . — Yield regions determined by boundary conditions 
in plane stress problems for circular plates. 
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Figure 3. — Yield surfaces for von-Mises and Hay thornthwai te yield 
conditions for plane stress problems. 




Figure 4. - Yield surfaces for Tresca and Johansen yield conditions 
for plane stress problems. 
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Figure 5* — General notion of a circle of di scontinuty . 




Figure 7. - General tine variation allowed for example 


eases. 
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1.0 



1.2 


96 



Figure 10. - Velocity distributions for various times for example 




Figure 12. — Velocity distributions for various tines for example 







